Critical behavior and entanglement of the random transverse-field Ising model 

between one and two dimensions 
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We consider disordered ladders of the transverse-field Ising model and study their critical prop- 
erties and entanglement entropy for varying width, w < 20, by numerical application of the strong 
disorder renormalization group method. We demonstrate that the critical properties of the ladders 
for any finite w are controlled by the infinite disorder fixed point of the random chain and the 
correction to scaling exponents contain information about the two-dimensional model. We calculate 
sample dependent pseudo-critical points and study the shift of the mean values as well as scal- 
ing of the width of the distributions and show that both are characterized by the same exponent, 
i'{2d). We also study scaling of the critical magnetization, investigate critical dynamical scaling 
as well as the behavior of the critical entanglement entropy. Analyzing the lo-dependence of the 
results we have obtained accurate estimates for the critical exponents of the two-dimensional model; 
p{2d) = 1.25(3), x{2d) = 0.996(10) and tp{2d) = 0.51(2). 



I. INTRODUCTION 

In nature there are materials, which are in a way be- 
tween two integer dimensions, such as they are built from 
(d — l)-dimensional layers having a finite width, w. Ex- 
amples are thin filmsi, magnetic multilayer^ or ladders 
of quantum spins^. One interesting question for such 
multilayer systems is the properties of critical fluctua- 
tions, when the linear extent of the layers, L, goes to 
infinity. If the system is classical having thermal fluctu- 
ations, finite-size scaling theorjiii^ can be applied. One 
basic observation of this theory is that for any finite w 
the critical behavior is controlled by the fixed point of 
the (d— l)-dimensional system, but the scaling functions 
in terms of the variable, w/L, involve also the critical 
exponents of the d-dimensional system. For example the 
critical points, Tc{w), measured at a finite width, w, ap- 
proach the true d-dimensional critical point, Tc = Tc (oo), 



Tc - Tc{w) 



w 



(1) 



where the shift exponent, Vg, generally corresponds to 
the correlation-length exponent, v, in the d-dimensional 
system. 

In a quantum system having a quantum critical point 
at zero temperature, T = 0, by varying a control pa- 
rameter, 9, the dimensional cross-over is a more subtle 
problem. If the d-dimensional critical quantum system is 
isomorphic with a, {d + l)-dimensional classical system^, 
then results of finite-size scaling can be transferred to the 
quantum system, too. This is the case, e.g. for the quan- 
tum critical point of the d-dimensional transverse-field 
Ising model, which is equivalent to the critical point of 
the classical {d + l)-dimensional Ising model. However, 
the situation is more complicated for antiferromagnetic 
models with continuous symmetry, such as for Heisenberg 



antiferromagnetic spin ladders. In this case the form of 
low-energy excitations could sensitively depend on the 
value of w: if the ladder contains even number of legs 
there is a gap, whereas for odd number of legs the sys- 
tem is gapless'^.. In the following for quantum systems we 
restrict ourselves to models with a discrete symmetry, 
such as to the transverse-field Ising model. 

In disordered systems, in which besides deterministic 
(thermal or quantum) fiuctuations there are also disor- 
der fluctuations in a sample of finite width one can define 
and measure a sample-dependent pseudo-critical point, 
Tc{w) (or 6c{w)), and study its distribution^. In particu- 
lar one concerns the shift of the mean value, Tc(u'), and 
the scaling of the width of the distribution, ATc{w). In 
this case besides the shift exponent, Vg, which is defined 
analogously to Eq. ((T|) one should determine the width ex- 
ponent, Viu, too, which is deflned by the scaling relation: 



(2) 



According to renormalization group theory^ the finite- 
size scaling behavior of random classical systems depends 
on the relevance or irrelevance of the disorder—. If the 
disorder represents an irrelevant perturbation at the pure 
system's fixed point, which happens if the correlation 
length exponent in the pure system satisfies Vp > 2/d, 
than for the disordered system we have Vs — and — 
2/d and the thermodynamic quantities at the fixed point 
are self-averaging. On the contrary for relevant disorder, 
which happens for Vp < 2/d, there is a new conventional 
random fixed point with a correlation- length exponent, 
V > 2/d^'^ , and we have Vs = = v- In this fixed point 
there is a lack of self-averaging. These predictions, which 
have been debated for some timeii, were checked later for 
various models^'^ i^^'^'^i^'* . 

For quantum systems quenched disorder is perfectly 
correlated in the (imaginary) time direction, therefore 
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generally it has a more profound effect at a quantum crit- 
ical pointi^. In some cases the critical properties of the 
random model are controlled by a so called infinite disor- 
der fixed point "'^^j in which the disorder fluctuations play 
a completely dominant role over quantum fluctuations. 
This happens, among others for the random transverse- 
field Ising model, as shown by analytical results^^ in Id 
and numerical result o^^i^^'^^i^° in 2d. Finite-size scaling 
has been tested for the Id model and a new scenario is 
observed^i. The finite-size transition points, denoted by 
Oc{L) in a system of length, L, are shown to be charac- 
terized by two different exponents, < Py^,. This means, 
that asymptotically A9c{L)/[9c{oo) — 9c{L)] — > oo, which 
is just the opposite limit as known for irrelevant disorder. 

In the present paper we go to the two-dimensional 
problem and study the finite-size scaling properties of 
ladders of random transverse-field Ising models. For 
this investigations we use a numerical implementation 
of the so called strong disorder renormalization group 
method^. As in 2d this method is expected to be asymp- 
totically exact in large scales. In the numerical imple- 
mentation of the method we have used efficient com- 
puter algorithms and in this way we could treat ladders 
with a large number of sites: we went up to lengths 
L = 4096 for w = 20 legs and used 4 x 10^ random 
samples. Our aim with these investigations is threefold. 
First, we want to clarify the form of finite-size scaling 
valid for this random quantum model. Second, using 
the appropriate form of the scaling Ansatz we want to 
calculate estimates for the critical exponents of the 2d 
model. Previous studiesiiii2iii^'22ii2^ in this respect have 
quite large error bars and we want to increase the ac- 
curacy of the estimates considerably. Our third aim is 
to calculate also the entanglement entropy^"' in the lad- 
der geometry and study its cross-over behavior between 
one^^ and two dimensions^SiiSl. 

The structure of the rest of the paper is the follow- 
ing. The model and the method of the calculation is pre- 
sented in Sec. [ITl In Sec. IIIII finite-size transition points 
are calculated and their distribution (shift and width) 
is analyzed. In Sec. IIVI we present calculations at the 
critical point about the magnetization and the dynam- 
ical scaling behavior. Results about the entanglement 
entropy are presented in Sec|Vl Our paper is closed by a 
discussion. 



II. MODEL AND METHOD 

A. Random transverse-field Ising ladder 

We consider the random transverse-field Ising model 
in a ladder geometry in which the sites, i and j, are 
taken from a strip of the square lattice of length, L, and 
width, w. We use periodic boundary conditions in both 



directions. The model is defined by the Hamiltonian: 

W = -^J,,afa^--^/»,af (3) 

in terms of the Pauli-matrices, cr^'^. Here the first sum 
runs over nearest neighbor sites and the Jij couplings and 
the hi transverse fields are independent random numbers, 
which are taken from the distributions, p{J) and q{h), 
respectively. For concreteness we use box-like distribu- 
tions: p{J) = 1, for < J < 1 and p{J) = 0, for J > 1; 
q{h) = 1/ho, ior < h < Iiq and q(h) — 0, for h > hg. 
We consider the system at T = and use 6 — In ho as 
the quantum control parameter. 

In the thermodynamic limit, L oo, the system in 
Eq.Q displays a paramagnetic phase, for 9 > 9c{w), and 
a ferromagnetic phase, for 9 < 9c{w). In between there 
is a random quantum critical point a.t 9 = 9c(w) and we 
are going to study its properties for various widths, w. 



B. Strong disorder renormalization group method 

The model is studied by the strong disorder renor- 
malization group method22,, which has been introduced 
by Ma, Dasgupta and Hu2^ and later developed by D. 
Fisher— and others. In this method the largest local 
term in the Hamiltonian (either a coupling or a trans- 
verse field) is successively eliminated and at the same 
time new terms are generated between remaining sites. 
If the largest term is a coupling, say J2.3 = O connect- 
ing sites 2 and 3, (fi being the energy-scale at the given 
RG step), then after renormalization the two sites form 
a spin cluster with an effective moment 112.3 — fJ-2 + fJ-3, 
where in the starting situation each spin has unit mo- 
ment, Hi = I. The spin cluster is put in an effective 
transverse field of strength: /i2,3 ~ ^2^3/>/2,3, which is 
obtained in second order perturbation calculation. On 
the other hand, if the largest local term is a transverse- 
field, say /i2 = f^, then site 2 is eliminated and new cou- 
plings are generated between each pairs of spins, which 
are nearest neighbors to 2. If say k and / are nearest 
neighbor spins to 2, than the new coupling connecting 
them is given by: Jkj ~ J2,kJ2,i/h2, also in second order 
perturbation calculation. If the sites k and I are already 
connected by a coupling, Jk.i ^ 0, than for the renormal- 
ized coupling we take max[Jk^i, Jk,i]- This last step is 
justified if the renormalized couplings have a very broad 
distribution, which is indeed the case at infinite disorder 
fixed points. The renormalization is repeated: at each 
step one more site is eliminated and the energy scale is 
continuously lowered. For a finite system the renormal- 
ization is stopped at the last site, where we keep the 
energy-scale, fl*, and the total moment, fi*, as well as 
the structure of the clusters. 
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C. Known exact results in the chain geometry 

The renormahzation has special characters in the chain 
geometry, i.e. with w = 1. In this case the topol- 
ogy of the system stays invariant under renormahzation 
and the couplings and the transverse fields are dual vari- 
ables. From this follows that at the quantum critical 
point the couplings and the transverse fields are deci- 
mated symmetrically, thus the critical point is located at 
6'c(l) = 0^. The RG equations for the distribution func- 
tion of the couplings and that of the transverse fields can 
be written in closed form as an integro-differential equa- 
tion, which has been solved analytically both at the crit- 
ical point^^ and in the off-critical region, in the so-called 
Griffiths-phase^. Here we list the main results. 

The energy-scale, fi, and the length-scale, L, are re- 
lated as: 

ln(r!o/0) ^ , (4) 

with an exponent: ipiXd) = 1/2. (Here L can be the 
size of a finite system and JIq is a reference energy scale.) 
The average spin-spin correlation function is defined as 
G{r) = [(erf CTf_|_^)]av, where (...) denotes the ground- 
state average and [. . . ]av stands for the averaging over 
quenched disorder. In the vicinity of the critical point 
G{r) has an exponential decay: 

G(r) ~ exp(-r/0 , (5) 

in which the correlation length, ^, is divergent at the 
critical point as: 

e~l0-0er^ (6) 

with v{\d) = 2. At the critical point the average corre- 
lations have a power-law decay: 

G{r)^r~^\ 6 = 9,, (7) 

with a decay exponent: 

x{ld) = {?>- \/5)/4 . (8) 

The average cluster moment, /j,, is related to the energy- 
scale, as: 

^i ~ M^o/m^ , (9) 

with (j>{ld) = (1 + %/5)/2. The average cluster moment 
can be expressed also with the size as /i ~ L'^f , where 
the fractal dimension of the cluster is expressed by the 
other exponents as: 

df = (P^p = d- X , (10) 

with d = 1. 



III. FINITE-SIZE CRITICAL POINTS 
A. Results in the chain geometry 

In the chain geometry finite-size critical points are 
studied in Ref.^^, in which they are located by different 
methods, which all are based on the free-fermion map- 
ping of the probleni^. The finite-size critical points are 
shown to satisfy the micro-canonical condition: 

L L 

^InJ, -^In/i, , (11) 

1=1 1=1 

from which follows that the distribution of 9c{L) is 
Gaussian with zero mean and with a mean deviation of 
A9c{L) ~ i"^/^. Consequently the width-exponent of 
the distribution is given by: 

= u{ld) = 2 . (12) 

On the other hand the shift-exponent is given by: 

i^s^l, (13) 

although in some cases (c.f. for periodic boundary con- 
ditions) the prefactor of the scaling function can be van- 
ishing. 

B. The doubling method 

In the ladder geometry, i.e. for w > 2, the free- 
fermionic mapping is no longer valid, therefore new meth- 
ods have to be utihzed to locate pseudo-critical points. 
Here we used the doubling method combined with the 
strong disorder renormahzation group. 

In the doubling procedureJ'^ for a given random sam- 
ple (a) of length L and width w, we construct a repli- 
cated sample (2a) of length 2L and width w by gluing 
two copies of (a) together and study the ratio of the 
magnetizations: rm(a, L, w) = m(2a, 2L, w)/m(a, L, w), 
which are calculated by the strong disorder renormahza- 
tion group method. In Fig[T] we illustrate the 9 depen- 
dence of the total magnetic moments, /i(2a, 2L, w) and 
fi{a,L,w), for a given sample of a w = 2-leg L = 128 
ladder. The corresponding ratio of the magnetizations, 
rm{a, L, w), is shown in the upper inset of this figure. It 
is seen, that in the ordered phase: 9 < 9c{a,L,w) this 
ratio approaches rm{a, L,w) 1. On the other hand 
in the disordered phase: 9 > 9c{a, L,w) the magnetiza- 
tions approach their minimal values, which in the SDRG 
method can be 1/2L and 1/L, respectively, thus we have 
rm(a, L,w) — *■ 1/2. In between there is a sudden change 
in the value of this ratio, which can be used to define a 
sample-dependent pseudo-critical point, 9c{a, L,w). 

There is another possibility, if we consider the ratio 
of the two gaps: rfi{a,L,w) = Q{2a,2L,w)/Q{a, L,w), 
which are also calculated by the strong disorder renor- 
mahzation group method. In FiglT] we show the two 
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log-gaps, — logri(2a, 2L, w) and —log^l{a,L,w), for the 
same sample as before and the corresponding ratio, 
rn{(x, L,w), is put in the upper inset of this figure. It is 
seen that this ratio in the ordered phase, 9 < 9c{oi, L, w), 
approaches rfi{a, L,w) and in the disordered phase: 
9 > 9c{oL, L,w), goes to ra{a,L,w) 1. In between 
this ratio has a quick variation and we can fix the point 
where ro(a,L,w) = 1/2 to define a sample-dependent 
pseudo-critical point, 9c{a, L, w). 



C. Numerical results 

1. Comparison of the two definitions 

In the doubling method we have calculated pseudo- 
critical points by using both ratios. Wc have observed, 
that for a given sample 9i^^ calculated from the ratio of 
the gaps is always somewhat smaller, than 9c"^\ which 
is obtained from the ratio of the magnetizations. This 
is illustrated in the upper inset of Fig[T]for a given sam- 
ple. We have also calculated for several realizations the 
ratio of the two pseudo-critical points, 9^"^^ /9^"^\ which 
is shown in the lower inset of Fig[T] as a function of 
for the w = 10 leg ladder for various lengths, L = 32, 64 
and 128. The relative difference between the two pseudo- 
critical points is indeed vary small, it is of the order of 
10"'^ and this is decreasing with increasing L and w. In 
the following we restrict ourselves to those pseudo-critical 
points, which are calculated from the ratio of the magne- 
tization and which have a relative precision of 10^** for 
each sample. 



2. Distribution of finite-size critical points 

We have calculated the distribution of pseudo-critical 
points for ladders with a fixed number of legs, 1 < w < 
20, for varying lengths, L ~ 2\ with / = 5, 6, . . . , 10. In- 
deed, for the largest values of L the relation, w/L <^ 1 is 
well satisfied. For the w = 10 leg ladder the distribution 
of the 9c values for various lengths are shown in Figl^l 
which are obtained for 10* realizations for each cases. As 
seen in this figure the width of the distribution is decreas- 
ing with increasing L and there is only a weak shift of 
the position of the maximum. The distributions some- 
what deviate from Gaussians, they are asymmetric, as 
can be seen in the log-lin plot in the inset of FigO With 
increasing L, however, the skewness of the distribution is 
decreasing, which is in agreement with the expectation, 
that in the w/ L — > limit we get back the corresponding 
results for chains. 
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FIG. 1: (Color online) SDRG results for the total mag- 
netic moments, /i(2a) and /i(a), as well as for the log-gaps, 
— logn(2a) and — logf2(a), as a function of the control pa- 
rameter, 9, for a given realization (a) of a, w = 2-leg L = 128 
ladder and its double (2a), see text. Upper inset: ratio of the 
magnetizations, r™, and that of the gaps, rn, as a function 
of 8 (log-lin scale) in the vicinity of the finite-size transition 
points. The finite-size critical point for is given at the 
jump, for rn it is located where rn = 1/2. Lower inset: Ratio 
of the two pseudo-critical points, rj = S^^'/Sc'"' as a function 
of si™' for the w = 10 leg ladder for various lengths and for 
500 realizations. 




FIG. 2: (Color online) Distribution of the pseudo-critical 
points, 8c, for the w — 10 leg ladder for various lengths and 
for 10^ realizations. In the inset in a log-lin plot deviations 
from the Gaussian distributions are seen, which however are 
decreasing with increasing L. 



3. "True" critical points for ladders 

For a fixed value of the number of legs, w, we have 
calculated the mean value of the pseudo-critical points. 
We have observed that the L-dependence of 9c{w, L) be- 
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TABLE I; Quantum critical points of ladders of the random 
transverse- field Ising model, 9c, and the asymptotic prefactor 
of the standard deviation in Eq.((T7}, a, for different number 
of legs. 



w 


6c 


a 


w 


6c 


a 


1 


0.00021(30) 


1.413(8) 


11 


1.39399(20) 


0.728(5) 


2 


0.64418(15) 


0.997(1) 


12 


1.41211(10) 


0.709(3) 


3 


0.94736(10) 


0.925(5) 


13 


1.42778(10) 


0.696(2) 


4 


1.08059(15) 


0.881(2) 


14 


1.44165(20) 


0.681(4) 


5 


1.16859(10) 


0.844(5) 


15 


1.45397(10) 


0.669(4) 


6 


1.23207(15) 


0.815(5) 


16 


1.46472(15) 


0.658(4) 


7 


1.27962(15) 


0.807(4) 


17 


1.47445(10) 


0.650(4) 


8 


1.31727(20) 


0.781(3) 


18 


1.48332(10) 


0.640(1) 


9 


1.34787(10) 


0.763(3) 


19 


1.49095(30) 


0.626(4) 


10 


1.37270(15) 


0.733(2) 


20 


1.49855(15) 


0.620(2) 



tern, 9c{2d), see FiglS) Here the corrections for large-w 
are expected to have a power-law form, and analogously 
to Eq.(IT]), it contains the shift exponent, i^s, of the 2d 
system. 

Estimates for the effective (w-dependent) values of the 
shift exponent are obtained from the ratio of the second 
and the first finite differences: 



1 



1 



(14) 



which are calculated at the central point of five-point 
fits. The effective exponents are given in the upper inset 
of FiglSl which are extrapolated as l/i'si'^d) = 0.81(10), 
thus we obtain: 



i/s(2d) = 1.24(15) . 



(15) 



comes weaker and weaker with increasing L, which is 
in agreement with the fact, that the system approaches 
more and more the chain geometry. Due to this one can 
obtain accurate estimates in the thermodynamic limit for 
the "true" critical points of ladders, which are listed in 
Table |T] for different number of legs. Here the errors are 
merely due to disorder fluctuations since for L > 192 the 
finite length effects are negligible. 

Here we also list our estimate for the chain, w — I, 
which agrees within the error of the calculation with the 
exact resuh: 61^(1) = and a(l) = \/2 (see Sec lHIC4l) . 



a> 1 




15 20 



FIG. 3: (Color online) Critical points of the ladders for vary- 
ing number of legs, w, and the extrapolation curve (broken 
(green) line) for large w. Upper inset: estimates of the in- 
verse of the local shift exponent, l/iya{w), calculated through 
Ea. (ll4|l . The broken (green) straight line indicates the ex- 
trapolation through 1/w. Lower inset: ratio of the scaled 
difference of the critical points and the scaled standard devi- 
ations as a function of l/w(see text). The horizontal broken 
(green) line at rs = 1. indicates a value which is close to the 
expected asymptotic behavior. 

These data approach the critical point in the 2d sys- 



4-. Scaling of the width of the distribution 

We have measured the standard deviation of the distri- 
bution of the pseudo-critical points, A6'c(w,i), for lad- 
ders with w legs and with a varying length, L. This 
quantity is expected to scale with the length as: 



(16) 



where the scaling function, cr(y), behaves for small argu- 
ments as: a{y) - y-i/>^u,(2d)+i/u,,{id) ^ p^.^^ ^^lis follows, 

that for finite widths we have: 



A0c{w,L) = L 



a{w) 



(17) 



with a prefactor, a(w), which behaves for large w as 
a(w) ~ w'^, with an exponent e — —1/1^^(2^) -I- l/^'u,(ld). 
We have checked this scenario by analyzing the data for 
A0c(w, L). First, for a fixed w we have fitted a func- 
tion a{w)L^^ , with a free parameter, u. We have found 
that for each widths, 1 < w < 20, the exponent oj agrees 
with \/vw{^d) = 0.5, within a few thousands of error, as 
illustrated in FigH) 

In the next step we have fixed the value of w — 0.5 
and estimated the limiting value of A9c{w, L)L'^-^ for 
large L, which is denoted by a(w). These limiting val- 
ues are presented in Table [H which are analyzed for large 
w. As seen in the upper inset of Fig. 0] in a log-log 
plot the a(w) values are asymptotically on a straight 
line. We have calculated effective, w-dependent expo- 
nents: e{w) = log(a(w)/a(u'/2))/log2, which are pre- 
sented in the lower inset of FigH] as a function of 1/w. 
These effective exponents have a weak w-dependence and 
we estimate its limiting value as e = —0.30(2). With this 
we have for the width exponent in 2c?: 



uU2d) = 1.25(3) 



(18) 



Closing this section we try to decide in a direct way 
about the relation between the two exponents, i^s(2d) and 
ly^i'^d). For this we form the scaled difference: d9c{w) = 
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FIG. 4: (Color online) Scaling of the width of the distribution 
of the pseudo-critical points, A9c{w, L), with L for different 
number of legs, w. In a log- log plot the slope of the curves 
are compatible with the theoretical prediction, Vw{ld) — 1/2, 
which is indicated by a full straight line. In the upper inset 
the limiting value of the prefactor, a{w) is shown as a function 
of w in a log-log plot. The dashed (green) straight line has 
a slope, e — 0.30, as extrapolated from effective exponents in 
the lower inset. 



Ai9c{w)w (see Ea.([T4|). which scales as w^^^'^", as well 
as the scaled standard deviation: sa{w) — w~^-^a{w), 
which scales as w~^^'^™, and form their ratio, rs(w) = 
dO c{w) / sa{w) . As seen in the lower inset of Fig [3] this 
ratio approaches a finite value which can be estimated 
as rs — 1.01(2). Thus we can conclude that at the infi- 
nite disorder fixed point of the 2d random transverse- field 
Ising model the shift and the width exponents are equal 
and they correspond to the correlation length exponent 
of the model. 

Using the best estimate for Vs{2d) — Vw{2d) and 
including the first analytic correction to scaling term: 
= edw) - ylw"i/'"=(l + B/w) we fit our data (see 
Figl3]) and obtain for the critical point of the 2d system: 

ec{2d) = 1.676(5) . (19) 

This value is in agreement with the previous estimate, 
ec[2d) = 1.680(5), in Ref.[13. 



IV. SCALING AT THE CRITICAL POINT 

Having estimates for the critical points of random lad- 
ders with w legs, 9c{w), we have calculated scaling of the 
magnetization at the critical point as well as the criti- 
cal dynamical scaling. These calculations are made for 
lengths up to 2^^ and for 4 x 10^ realizations. 



A. Magnetization 

We have calculated the average total magnetic moment 
at the critical point, ^c{w,L), for varying lengths, L, 
which is expected to scale as: 

fi,{w,L)=L^f^^^^fi,{w/L) , (20) 

with a scaling function, which behaves for small argu- 
ments as: iic{y) ~ yd}(2d)~df{id) ^ Then, for a finite 
width, w, we have: 

^i,{w,L) = L^i^^^^h{w) , (21) 

with a prefactor, which for large w behaves as: b{w) ^ 
w'^, with K = df{2d) — d/(ld). 
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FIG. 5: (Color online) Scaling of the average total magnetic 
moment at the critical point of a ladder with w legs and length 
L, ^c{w,L). In the log-log plot the slope of the curves are 
compatible with the theoretical prediction, d/(ld) in Eq.((TO}, 
which is indicated by a full straight line. In the lower inset 
the limiting value of the prefactor, b{w) is shown as a function 
of w in a log- log plot. The dashed (green) straight line has a 
slope, K — 0.195, as extrapolated from effective exponents in 
the upper inset. 

The scaling Ansatz in Eq. ([2T|) is checked in FiglS] 
Then, we have calculated the limiting value of 
L~'^f^^'^^ lj,c{w, L), which is denoted by b{w) and which 
is presented as a function of w in a log-log plot in the 
lower inset of Figl5l Effective, w-dependent exponents 
are calculated, which are extrapolated in the upper inset 
of Fig O giving k = 0.195(10). Thus the fractal dimension 
in 2d is df{2d) ~ df{\d) + n and from Eq.ijTU]) we obtain 
for the magnetization scaling dimension: 

x{2d) = 0.996(10) . (22) 

B. Dynamical scaling 

At an infinite disorder fixed point there is a special 
form of dynamical scaling, as given in Eq.(|3]). The en- 
ergy scale of a sample at the end of the renormalization 
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then: 
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FIG. 6: (Color online) Distribution of the last decimated 
log-coupling, A, (upper panel) and the last decimated log- 
transverse field, 7 (lower panel) for various lengths, L, for the 
w = 10-leg ladder. In the insets the distribution of the scaled 
variables: XL''^'^ and -fL'"'^, respectively are shown. 



T{w,L) ^ L'^^^'^'>g{w) , 



(24) 



with glw) ~ for large w with 6 — ^p{2d) — ip{ld). 
2.5 
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FIG. 7: (Color online) Scaling functions for typical energy- 
scales: g{w) = lim_L_cx) G(w, L)L~^'^^, in which G{w,L) is 
either T{w,L), Ar(iu,L), A{w,L) or AA{w,L), see text. 



We have checked that the scaling form in Eq. ([Ml) is 
indeed satisfied for all values of 1 < w < 20 and than 
calculated the limiting value of r(w, which is 

denoted by g{w). As illustrated in Figl7]the scaling func- 
tions of the typical energy-scales have only a very weak 
w dependence, and we estimate (not shown) a small ex- 
ponent: S — 0.01(2). Thus we have for the 4' exponent 
in the 2d model: 



tp{2d) = 0.51(2) 



(25) 



can be defined either by the value of the last decimated 
(log) coupling — log J = A or by the last decimated (log) 
transverse field — log /i = 7. The distribution of A as well 
as 7 are shown in FiglHlin upper and in the lower panel, 
respectively, for the w — 10-leg ladder for various values 
of the length, L. An appropriate scaling collapse of the 
date is observed in terms of the scaling variables, \L~^ 
and 'yL~^, with ip = tp{ld) = 1/2, as illustrated in the 
insets. 

In order to have a more quantitative picture about dy- 
namical scaling we consider the mean value: T{w, L) = 
[j{w, L)]av and the standard deviation, Ar(ii;, L) and 
similarly, A{w,L) = [A(w,L)]av and AA(ii;, L). All these 
quantities are expected to scale in the same way, for ex- 
ample with T{w,L) we have: 

r(w,L) =i'^(2d)p(^/^) ^ (23) 
with f(y) - y^{-2d)-4,{id) ^ Yov a finite width, w, we have 



V. ENTANGLEMENT ENTROPY 

In the ladder geometry we consider a block. A, which 
contains all the w legs and has a length, £ L. Con- 
sequently the block has two parallel lines of width, w, 
at which it has contact with the rest of the system, B. 
The entanglement of A with B is quantified by the von- 
Neumann entropy^!: 

5^(u.,^) = -Tr^(p^logp^) , (26) 

in terms of the reduced density matrix: p_4 = Trgj^') {^\, 
where I^E*) denotes a pure state (in our case the ground 
state) of the complete system. 

At the critical point of a random quantum system the 
properties of which are controlled by an infinite disor- 
der fixed point the asymptotic behavior of the entropy 
in the large I limit can be obtained by the strong dis- 
order renormalization group method. For the random 
transverse-field Ising model entanglement between A and 
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B are given by such renormalized spin clusters, which 
contain sites both in A and in B, and the cluster is elim- 
inated at some point of the renormalization^^^i^l. Due 
to the very broad distribution of the effective couplings 
and transverse fields, the cluster at the energy scale of 



its decimation is in a so called GHZ entangled state of 
the form: l/v^d TT ■ ■ • T) + I ii • ■ ■ i))- Each such clus- 
ter contribute by an amount of log 2 to the entanglement 
entropy, thus calculation of the entropy is equivalent to 
a cluster counting problem, which is illustrated in FiglH 



J 




FIG. 8: (Color online) Structure of the decimated spin clusters at the critical point of a ladder of 20 x 128 sites, which is 
devided into two equal blocks the boundary of which is indicated by thick vertical lines. There are three clusters, denoted by 
blue, red and green colors, respectively, which contain sites at both blocks and thus result in an entropy 3 x log 2. 



In the chain geometry the asymptotic behavior of the en- 
tropy at the critical point is obtained from the analytical 
solution of the RG equations as^: 

SAl,^)-^iog£ + k{l) , (27) 

where k(l) is a non-universal constant, which depends 
on the form of the disorder, whereas the prefactor of 
the logarithm, c(l), which is also called as the effec- 
tive central charge, is universal and given by: c(l ) = 
In 2/2. This result is checked numerically in Ref . [33| . In 
the two-dimensional case, which is expected to hold for 
w/L — 0(1), there are somewhat conflicting numerical 
results at the critical point. Lin et al^ have observed a 
double-logarithmic multiplicative factor to the area-law: 
'Sa{^i^) ~ ^loglog.^ whereas later Yu et al.^'^ argued to 
have only a subleading logarithmic term to the area law: 
5^(£,£) a£ + b\ogi + k. 

Here we study numerically the critical ladder systems 
with various number of legs and try to identify the cross- 
over between one- and two dimensions. To illustrate the 
I dependence of the entanglement entropy in Fig[9] we 
show 5^ as a function of log^ for different number of 
legs for L = 4096. (We have checked, that the asymptotic 
results does not change for L = 2048.) The central parts 
of the curves are very well linear having approximately 
the same slope, which is consistent with the exact result 
for the w ^ I chain geometry. Thus we conclude that 
the effective central charge, c(w), does not depend on 
the number of legs. 

In the next step we fix c(w) = In 2/2, calculate the non- 
universal term: k{w, £) = Sy\^{w, i) — ^ log£ and take its 
limit, k{w), for large i (but still with £ <C i). As illus- 
trated in the upper inset of Fig[5]the ^-dependent correc- 
tion term is approximated as log^/£ and the asymptotic 




10"'' 10° 10\ 10^ 10^ 



FIG. 9; (Color online) The log-^-dependence of the entan- 
glement entropy at the critical point of ladders for different 
number of legs, w, and for L = 4096. The linear part of the 
curves has approximately the same slope, which is consistent 
with In 2/6, as indicated by a full straight line. Upper in- 
set: the non-universal part of the critical entropy, k(w,t), for 
the w = 10-leg ladder and its extrapolation for £ S> if (but 
£ <^ L — 4096) with a correction term of ~ log£/£. Lower 
inset: the asymptotic value of the non-universal part of the 
critical entropy, k{w), as a function of w. Asymptotically 
there is a linear tii-dependence, which is shown by a broken 
(green) straight line. 



non-universal terms, k{w), are shown for different num- 
ber of legs in the lower inset of Fig[9l One can see that 
starting with the chain, w = 1, first k{w) is decreasing, 
has a minimum around w — 3 and then starts to increase. 
This increase for large w is approximately linear, we have 
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fitted: k{w) = 0.0256(5)u;+0.148(5). This linear increase 
is compatible with the area law, which should hold for 
non-critical systems and for large blocks. Our analysis 
can be used to clarify the one- to two-dimensional cross- 
over of the entropy in the limit w/£ ^ 1. However, our 
data can not be used to make predictions further, for 
w/L = 0(1), i.e. for the two-dimensional case. For this 
one should analyze the occurrent and possibly very weak 
w dependence of the prefactor of the linear term of k{w), 
which however can not be done with our data, which are 
only up to w = 20. 



VI. DISCUSSION 

In this paper we have studied the critical properties 
and the entanglement entropy of random transverse-field 
Ising models in the ladder geometry by the strong dis- 
order renormalization group method. In our numerical 
calculation we went up to w=20 legs and with a length up 
to L=4096 for 4 x 10'' realizations. In principle the sizes 
of the systems could have been increased further, but it 
was not necessary. With L we have already reached the 
limit where no further systematic finite-size effects are 
seen. On the other hand for larger values of w we would 
have obtained too large errors in calculating quantities, 
such as through two-point fit. 

First, we have calculated sample dependent finite-size 
critical points, which are obtained by the doubling pro- 
cedure and the strong disorder renormalization group 
method. We have analyzed the shift of the mean value 
of the transition points and the width of the distribution 
as a function of the number of legs, w, and estimated 
the exponents of the 2d model, i^s(2(i) and i/u,(2(i), re- 
spectively. These are found to be identical and given by 
the correlation-length exponent of the 2d model. Conse- 
quently the scaling behavior of the pseudo-critical points 
of the 2d random transverse-field Ising model is in the 
same form as for classical and conventional random fixed 
points^. In this respect there is a difference with the Id 
model2i, in which Vsi^d) < Vwi^d). For this latter model 
probably the free-fermionic character could be the reason 
for the different scaling properties. Our estimate for the 
correlation length exponent, iy{2d) = 1.25(3), is clearly 
larger than the possible limiting value of 2/d, which has 
been observed in the Id model and in some other random 
systemsii. 

Scaling at the critical point for different quantities is 
analyzed in a similar way, what we summarize here as 
follows. Let us consider a physical observable. A, which 
at the critical point has the mean value, A{w^L). This 
quantity scales with the critical exponent of the 2d model, 
a{2d), as: 

A{w,L) ^ L'^^-^'^^iw/L) (28) 
where the scaling function, A{y), for small arguments 



TABLE II: Numerical estimates of the critical exponents at 
the infinite disorder fixed point in 2d. MC: Monte Carlo 
simulation; SDRG: numerical strong disorder renormalization 
group; CP: Monte Carlo simulation of the 2d random contact 
process. The exponents, <j), denoted by an asterisk are calcu- 
lated from the scaling relation in Eq. (|10|) . 



i, 







X 


method 


0.4(1) 


2.5* 




1.0 


MOii 


0.42(6) 


2.5(4) 


1.07(15) 


1.0(1) 


SDRGi^ 


0.5 


2 




0.94 


SDRGi^ 


0.6 


1.7 


1.25 


0.97 


SDRG^ 


0.51(6) 


2.04(28)* 


1.20(15) 


0.96(2) 


CP?^ 


0.51(2) 


1.97(10)* 


1.25(3) 


0.996(10) 


this work 



behaves as: 

A{y) ^ y"(2<i)-a(lrf) (29) 

where Q;(ld) is the critical exponent in the Id model. 
Consequently for a finite w, but for L — )■ oo, we have 

A(«;,L) - L"(i'^)a(w) (30) 

with a(w) ~ w'^ and ld — Q;(2d) — a(l(i). In general we 
measure the scaling function a(vj) for different widths, 
estimate the exponent w and calculate the critical expo- 
nent in 2d as: q;(2o?) = a(lc?) -I- lo. Since the exponents 
in \d are exactly known and the correction term, w, is 
comparatively small we have obtained quite accurate ex- 
ponents in 2d. In the following we compare the estimates 
for the different critical exponents in the 2d infinite dis- 
order fixed point, which are listed in Table [Til 

Here besides different numerical strong disorder renor- 
malization group results there are also Monte Carlo simu- 
lations, both for the random transverse-field Ising model 
and for the random contact process. This latter model 
is expected to belong to the same universality classes, at 
least for strong enough disorder. It is seen in Table |TT] 
that our estimates fit to the trend of the previous results 
and generally have a somewhat smaller error. 

We have also studied the scaling behavior of the en- 
tanglement entropy in the ladder geometry. For a fixed 
width, w, the entropy is found to grow logarithmically 
with the length of the block, ^, and the prefactor is found 
independent of w. On the other hand the i independent 
term of the entropy is found to have a linear w depen- 
dence, at least for large enough w, which corresponds to 
the are law for this systems. 

The investigations presented in this work can be nat- 
urally continued for larger and larger widths and ap- 
proaching the case, wjL = 0(1), which corresponds to 
the two-dimensional model. However, with increasing 
w the numerical computation becomes more and more 
costly. The reason for this is the fact that the con- 
nected clusters in the strong disorder renormalization 
group method are typically of size w x w, which for large 
w becomes fully connected after decimating a small per- 
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cent of the transverse fields. The number of further renor- 
malization steps grows in a naive approach as , so that 
by this method one can not go further than L ~ 100 
or 200 in 2d. To treat larger 2d systems improved al- 
gorithms are necessary. Studies in this direction are in 
progress. 
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